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Abstract 

The decoration or iteration transformation was widely applied to solve exactly the magnetic spin models in 
one-dimensional and two-dimensional lattice. The motif of this letter is to extend the decoration transformation 
approach for models that describe interacting electron systems instead of spin magnetic systems, one illustrative 
model to be studied, will be the spinless fermion model on diamond chain. Using the decoration transformation, 
we are able to solve this model exactly. The phase diagram of this model was explored at zero temperature as well 
as the thermodynamics properties of the model for any particle density. The particular case when particle-hole 
symmetry is satisfied was also discussed. 
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1 Introduction 

Exactly solvable models in statistical physics and mathematical physics is one of the most challenging topics. 
Recently several exactly solvable models were studied in quasi-one dimensional classical-quantum models such as 
Ising-Heisenberg models [TJ [21 121 HI as wen as two-dimensional Ising models and Ising-Heisenberg models [71 [3], 
that can be mapped onto exactly solvable vertex models. Although the decoration transformation was introduced in 
the fifty decade [9l [10], in order to study the decorated spin models. Due to a successful application of this approach, 
recently this method was extended for a general case of decoration transformational], later extended by Strecka[12j 
even for classical-quantum models such as Ising-Heisenberg models. Recently another interesting application of 
decoration transformation was also investigated by Pereira et al. |13j where they considered a delocalized interstitial 
electrons on diamond-like chain and they also investigate the magnetocaloric effect in kinetically frustrated diamond 
chain [13] . meanwhile Strecka et aZ.|15] discussed the localized Ising spins and itinerant electrons in two-dimensional 
models, as well as two-dimensional spin-electron with coulomb repulsion[16 . 

On the other hand studies of strongly correlated electron systems are certainly one of the most important 
areas of condensed matter physics. A large number of papers have been concerned with the investigation of heavy- 
fermion behavior, magnetism of strongly correlated systems, high temperature superconductivity, or metal-insulator 
transitions and charge-ordering phenomena. A limiting case of these systems could be the spinless fermion model, 
where the spin orientation could be ignored or it can be understand as a fully polarized system. The one-dimensional 
case of this model already was investigated by Czart[T7] and Zhuravlev et al.[19\. Although spinless fermion models 
can be transformed onto XXZ models in a magnetic field parallel to the anisotropy[T71 Q15] by the use of the well 
known Jordan- Wigner transformation, this kind of mapping could be more involving for the Hubbard-like systems, 
where the spin orientation is considered. Therefore, the purpose here is how to use the decoration transformation 
approach for interacting electron systems, without mapping onto spin models. 

The outline of this letter is as follow: In sec. 2 we introduce the model to be studied. In sec. 3 we present 
the phase diagram at zero temperature. In sec. 4 is devoted the decoration transformation and transfer matrix 
approach, in order to obtain its exact solution. Whereas in sec. 5 is discussed the thermodynamics, particle density 
and correlation function of the model considered, and finally in sec. 6 we present our conclusions. 
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Figure 1 : Schematic representation for spinless fermion model on diamond chain, by vertical solid line we represent 
a hopping term and Coulomb repulsion between two sites, while by dashed line we represent only the Coulomb 
repulsion term. 

2 The model 

The model we consider in this letter is the spinless fermion model on diamond chain as displayed in fig. [TJ where by 
vertical solid line we mean the particle hopping term and coulomb repulsion term, while by dashed line we represent 
only the coulomb term repulsion. Therefore the Hamiltonian for spinless fermion model on diamond chain (fig. [1]), 
can be expressed by 

N 

H = J2 H i,i+i, (1) 
with N being the number of unit cell, whereas ffj^+i is given by 

Hi,i+i = - t (",",./",.- + al ida^ - n \ n a ,i + + ^(«c,i + «c,i+i)^ 

V 

+ V\n a ^n h ^ + — (n Cti + n c ,i+i) (n a ,i + , (2) 

where a ay i(a i a 4 ) are Fermi annihilation (creation) operators for spinless (or completely polarized) fermion respec- 
tively, with a — {a, b, c}, while n a j = a* a f a ai j being number operators. The Hamiltonian parameter t is the hopping 
term (solid lines), fi is the chemical potential, and the Coulomb repulsion term between fermions on neighboring 
sites (dashed lines) is represented V and V\. We also consider a periodic boundary condition. 

The spinless fermion model on diamond chain, has an important particle-hole symmetry [18J. For the purpose 
of this discussion, the following canonical transformation i — > a a ^ and a a ^ — > a) a i was considered. For the 

occupation number operator the transformation reads n a ,i — > a a ,ia} a i — l — n a ^, thus the transformed Hamiltonian 
becomes 



H' 



=* ( a i,i a fc,i + a l,i a a,i) - (V + Vi- fj) (n 0)j + n 6 ,j) - (V - ^)(n c .i + n c>i+ t) 
V 

+ Vin a4 n b ,i + — (n c s + n c . i+ i) (n aA + n M ) + (Vi + 2V - 3ju) . (3) 

The different signs of the hopping term on sites a and b can be compensated by a further canonical transformation 
as follows. By formally replacing say a\ , — > —a} a i {a a i — > —a a i ) for a sites, but keeping unchanged the operators 
on sites 6, or vice-versa. 

On the other hand, using the Jordan- Wigner transformation it is possible to map the spinless fermion diamond 
chain onto a Ising-Heisenberg diamond chain with non- uniform external magnetic field, as discussed in reference [4, 5|, 
the z-component of the spin-1/2 operator is related to number operator by n — a z + 1/2, while the creation 

(annihilation) operators are related through at = 2 J ^Ilfc<j a k) a J an< ^ a j = ^ (llfc<i a t ■ Therefore the 
models are related through the relations J H = V\, Ji = V, Hh = — (V + Vi)/2, Hj = fi — V and Jff A = —t/2. 

The Hamiltonian of spinless fermion model on diamond chain, has not been discussed yet anywhere . In this sense 
discussing this model, we could be opening several variants of fermion models and the most interesting models are 
the extended Hubbard-like models[20J with decorated interactions, where is taking into account the spin orientation, 
that is why it is interesting first to discuss the spinless fermion model. 
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3 The phase diagram 



Considering the state vectors, where hopping term t is acting at sites a and b. In each state there are only two 
possibilities labeled by or 1, it means an empty or occupied particle state respectively, this leads to the following 
state vector |W) = ci|0, 0} + C2 10, 1) + c%\l, 0) + c^\l, 1), where c, are the coefficients to be determined for each state. 
The states acting on Hamiltonian ((2j) at sites a and b for each elementary cell becomes, 



Hi,i+l\ 



0,0) = 


- - i n c,i 


0,1) = 


(V p 
V 2 2 


1,0) = 


YV p 
V 2 2 


1,1) = 





(n. 
(n. 



c,i+i)|0,0), 



n Cji+ i) 



|0, 1> — t[l,0), 
|1,0> — 1[0, 1>, 



+ n 0ti+1 ) + Vi-2p |1,1) 



(4) 
(5) 

(6) 
(7) 



After diagonalize the above 4x4 matrix we have 4 eigenvalues, these eigenvalues depends only on the number of 
operators n c j, along with the Hamiltonian parameters. The states (|4|) and ([7|) are already in their eigenstates, while 
the eigenvalues of eqs. ([5]) and (j6]) is given simply by — ^) (n c ,i + ti C) j+i) — p] ± t, whereas their respective 
eigenvectors are given by \v) s<a — A= (|0, 1) ± |1,0)), symmetric and anti-symmetric states respectively. Although 
there are 16 eigenvalues only four possible ground states energies were found, whose eigenvectors for different phases 
are given by 



JY 



|£0)=JJ|0,0)i® |0)i, 

i=l 
N 

i 5i >=n^(i i ' o >< + i o ' i ><)®i°>*' 
i=i 

N 

i^ 2 >=n^(i i ' o ^+i o ' i w®i 1 ^ 

i=l 
N 

|53>=JJ|l,l>i®|l) 4 



p = 0, 

P = l, 
p = 2, 
p = 3, 



(8) 
(9) 
(10) 

(11) 



where the states of type \a, b) <g> |c) corresponds to the particles states at sites a, b and c, respectively. 

The ground state phase diagram of V as a function of /i is displayed in fig. [5J for fixed values of t = 1 and 
Vi = V. The four states are limited as follow: 



| SO) : limited by /x -1, with p = 0, 
|S1) : between /i ^ —1 and V = /i with p = 1, 
| S2) : between 7 = p and V = (p - l)/2 with p = 2, 
| S3) : between 7 = (p - l)/2 and V = with p = 3. 



(12) 
(13) 
(14) 
(15) 



The first one is |S0), this corresponds just to a simple empty lattice particle (or fully-filled holes) on diamond 
chain with total density p = (yellow region). The second state is represented by |S1), where one particle is at 
site a or b of the elementary cell, then the particle density for this state is p = 1 (green region). Furthermore 
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Figure 2: Phase diagram of spinless fermion model on diamond chain as a function of /i and V, for fixed values of 
t = 1 and Vi = V. 

there is another state 152), this state corresponds to the configuration that one particle is fluctuating at sites a 
and b, whereas on site c of elementary cell there is another particle, then the total density is p = 2 (red region). 
Finally there is a fully-filled particle diamond chain state represented by \S3) with density p — 3, or it can also be 
understood as an empty lattice of hole state (blue region) . 

4 Decoration transformation and transfer matrix 

In order to study the thermodynamics properties of the spinless fermion model on diamond chain, it will be used 
the decoration transformationPJJ [TTJ [T^] as described bellow. Actually, it is not necessary to map the spinless 
fermion model onto spin models, like Ising-Heisenberg model[TJ [2j [3l IU [6]. Thus, there is an interest in performing 
the decoration transformation for operators. In this sense it will be apply directly the decoration transformation 
approachpjjl [Til 02] for the proposed model. The main aim to solve the Hamiltonian (jTJ) of spinless fermion model 
on diamond chain is to map onto an effective spinless fermion model without hopping term, whose Hamiltonian is 
given by 

N 

where V and fl are coefficients to be determined using decoration transformation[TO|, lllj. 
The Boitzmann factor of effective spinless fermion model can be expressed as follow 

w(n Cji , n C)i+ i) = cxp (^/3Vn Ctl n^ i+1 - f3^(n C:i + n C:i+ i) \ , (17) 

where /3 — 1/kT, with k being the Boitzmann constant and T the absolute temperature. 

On the other hand the Boitzmann factors for the spinless fermion model on diamond chain given for the 
Hamiltonian (JTJ) reads as 



Vn Cti n Cii+1 



2 K 



n 



c,i+\) 



(16) 
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w(n c ,i,n C;i+ i) = tr a , b (e^*^ 1 ) , (18) 
the operators n c> i ranges from to 1, then we have explicitly the Boltzmann factors, 

w(0, 0) =1 + 2c^ cosh(/3t) + e -P Vl+2 ^ 1 , (19) 
w(0, 1) =e^/ 2 + 2e^ 2 ^ v / 2 cosh(/3i) + e -P<Y+Vi)+*P*/2 t (20 ) 
«,(1, 1) = e ^ + 2e^"^ cosh(^) + e -^ 2V +^)+ 3 ^. (21) 

For simplicity the Boltzmann factor for diamond chain are denoted by wq = u>(0, 0), and w± = w(X, 0) = w(0, 1) 
and W2 — w(l, 1). 

In order to apply decoration transformation we need to impose the following condition: 

Z = fZ eff , (22) 

The Boltzmann factor for both systems must be equal in order to satisfy the decoration transformation ]!!)! 111] . 
Thus there are three unknown algebraic equations and three unknown parameters, this algebraic system is solved 
easily using the decoration transformation method, whose solutions are written as 

/ = «*, A=M-l V=lln(^_) (23) 



ft \WQ J P V 102 Wo 

where the /, ft and V are expressed as a functions of the original parameters of the Hamiltonian by means of wq, 
w\ and w 2 - 

On the other hand, the Boltzmann factor for effective spinless fermion model reads as 

w = 1, wi = e^ /2 = x, w 2 = e^-P 9 = x 2 y, (24) 

with being x = exp(/3ft/2) and y = exp(— /3V). 

In what follows, we are interested to solve the effective spinless fermion without hoping term (also known as the 
atomic limit), to solve this effective model it will be used the transfer matrix method[21J, given by 

T = ( I .fa V (25) 



x yx 

Note that the terms ft and V can be obtained in a similar way as were obtained for spin models. The eigenvalues 

2 -I- . I ( 1 — „,^ 2 j_ aZi 



of transfer matrix are given by A± = I 1 + yx 2 ± y (1 - yx 2 ) + Ax 2 J /2. 

Using the largest eigenvalues A+ of the transfer matrix T, we obtain the partition function per site of the model 
in terms of the effective spinless fermion model in atomic limit Z — f Z e ff, with Z e ff = A+ is the partition function 
for effective spinless fermion without hopping term. The partition function per elementary cell, is expressed in 
terms of the Boltzmann factors for the spinless fermion model on diamond chain, may be written 



Z = i (wo + w 2 + sj '(w - w 2 f + 4w 2 j . (26) 

From the partition function of spinless fermion model on diamond chain, it is possible to obtain the free energy 
by the relation fl = —kThiZ. Once known this result, we are ready to study several physical amounts such as 
entropy, specific heat, average energy and so on. 
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(a) V = 1 and Vi = 1 (b) V = 1 and Vi = 0.1 (c) V = 1 and V = 1 (d) V = 1 and Vi = 0.1 

Figure 3: In (a) and (b) is plotted the total density p as a function of chemical potential p, for fixed values of 
temperature. In (c) and (d) we display the amount (n a ) — (n c ) versus p, for same set of parameters. 



5 Thermodynamics and correlation functions 

So far, we have not yet specified the particles number of the diamond chain. In order to obtain the electron density 
per elementary cell, we take the derivative of free energy in relation to the chemical potential of spinless fermion 
on diamond chain 

= 2<n a ) + (n c ), (27) 



here we are assuming the exchange invariance of sites a and b. 

The particle density of type c can be obtained from effective spinless fermion model, by the relation 



P= dp 



1 



whereas the particle density for sites a and b per site, can be obtained combining eqs. (|2T|) and (|2"5)h resulting in 

(n a ) - P -P±. (29) 



Further we illustrate the chemical potential behavior for fixed values of t = 1, V = 1 and V\ = {1,0.1}. In fig. 
(3](a-b) we display the density p as a function of chemical potential p, for a range of temperatures T — 1.0, 0.4, 
0.1 and 0.01, this behavior is in agreement with phase transition at zero temperature (fig. ^ discussed previously. 
We show three plateaus at low temperature, when p becomes 1, 2 or 3. The size of plateau at density p = 2 is 
proportional the parameter V\, for small V% there is a short plateau (figj3jb)), while for large V\ long plateau is 
observed (the last one not illustrated). In fig. [HJc-d) the amount (n a ) — (n c ) versus p is displayed, in order to study 
the average particles number behavior at sites a and b compared to that one on site c. So for negative chemical 
potential we have (n a ) > (n c ), while for positive chemical potential we have (n a ) < (n c ). 

Another interesting properties we would like to discuss is the symmetry particle- hole described by eqs ([3]) must 
be satisfied only when V\ — V, therefore we find the following relation for the free energy, 



Cl(t, V, p) = Q(t, V, 2V - p) + (3V - 3p). (30) 

In general, the chemical potential depends on the temperature, however at half-filled band particle density 
p = 1.5 and under particle- hole symmetry relation, the chemical potential becomes independent of the temperature, 
given simply by p = V . If we look at the phase diagram displayed in fig. [5J the relation V = p corresponds the 
phase transition between density p = 1 and 2, as expected the density in this region should be p = 1.5, due to the 
thermal fluctuation begins to act. 
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Figure 4: For several values of density p is displayed the chemical potential as a function of temperature, for 
half-filled band. 

Furthermore, in figH] the chemical potential as a function of the temperature for fixed values of density p is 
displayed, from this plot we may conclude that there is a chemical potential independent on the temperature only 
at half-filled band and when symmetry particle-hole is satisfied. 

Another interesting quantity we would like to comment is the correlation function, it can also be obtained using 
transfer matrix method[21J, which after some algebraic manipulation we obtain, 



Some other nearest correlation or expected values of higher order can be obtained, combining the derivative of 
diamond chain Hamiltonian with respect to one of its parameter, and the decoration transformation for correlation 
function as discussed in reference [TUl E] . 

6 Conclusions 

Although the decoration transformation has originally been developed to study magnetic spin mo dels [IUJ El HU; an d 
since then, this method was widely applied to solve several decorated spin models. Motivated by this successful ap- 
plication, this transformation later was generalized [TJJ, and recently Strecka[12j extended even for Ising-Heisenberg 
spin models. However, there is no any approach developed for electron interacting systems, such as strongly corre- 
lated electron systems. In this sense, we discuss how the decoration transformation can also be used to map from 
one decorated electron system onto another effective electron system. As an illustrative application of decoration 
transformation for electron interacting models, it has been considered the spinless fermion model on diamond chain, 
where vertical solid line corresponds to hopping term and repulsive Coulomb interaction term, while by dashed line 
we mean the repulsive Coulomb interaction term between the nearest neighbor (fig. [lj. Furthermore, we discuss 
some properties of this model, such as the phase diagram at zero temperature, showing four different states with a 
given number of particles, changing from empty lattice (or fully-filled holes) of particles to fully-filled particles (or 
empty lattice of holes) on the diamond chain. The thermodynamics of this model, allows to display the particle 
density as a function of chemical potential per elementary cell as well as the chemical potential as a function of 
temperature. Finally an additional quantity also was considered such as the correlation function. 




(31) 
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